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ABSTRACT 


Experimental data from jet-engine tests have in- 
dicated that turbine efficiencies at takeoff can be as 
much as two points higher than those at cruise con- 
ditions. Recent studies have shown that Reynolds 
number effects contribute to the lower efficiencies at 
cruise conditions. In the current, study numerical sim- 
ulations have been performed to study the boundary 
layer development in a two-stage low-pressure tur- 
bine, and to evaluate the models available for low 
Reynolds number flows in turbomachinery. In a pre- 
vious study using the same geometry the predicted 
time-averaged boundary layer quantities showed ex- 
cellent agreement with the experimental data, but 
the predicted unsteady results showed only fair agree- 
ment with the experimental data. It was surmised 
that the blade count approximation used in the 
numerical simulations generated more unsteadiness 
than was observed in the experiments. In this study a 
more accurate blade approximation has been used to 
model the turbine, and the method of post- processing 
the boundary layer information has been modified to 
more closely r resemble the process used in the exper- 
iments. The predicted results show improved agree- 
ment with the unsteady experimental data. 


NOMENCLATURE 


C j Skin friction coefficient 

H Shape factor 
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M Mach number 

N Nozzle/vane 

P Pressure 

P t Total Pressure 

Re Reynolds number (axial chord, inlet vel.) 

S Entropy, arc length 

SSL Normalized suction surface arc length 

Tu Turbulence level 

u Local velocity 

b'oo Free stream velocity 

5* Displacement thickness/suction surface length 

rj Efficiency^ 

A Pressure gradient parameter 

i / Kinematic viscosity 

p Density 

cr Intermittency function 

9 Momentum thickness 

Q, Rotational velocity 

SUBSCRIPTS 


tt Total- to- total quantity 

wet Wetted distance 

9 Momentum thickness 

1 First nozzle inlet 

2 First nozzle exit/first rotor inlet 

3 First rotor exit/second nozzle inlet 

4 Second nozzle exit/second rotor inlet 

5 Second rotor exit 


INTRODUCTION 


Experimental data from jet-engine tests have in- 
dicated that unsteady blade row (wake) interactions 


and separation can have a significant impact on the 
efficiency of turbine stages. The effects of these in- 
teractions can be intensified in low-pressure turbine 
(LPT) stages because of the low Reynolds number 
operating environment. Measured turbine efficien- 
cies at takeoff can be as much as two points higher 
than those at cruise conditions [1]. Thus, during the 
last decade a significant amount of effort has been 
put into determining the effects of transit ion and tur- 
bulence on the performance of low pressure turbine 
stages. Experimental investigations have been per- 
formed, for example, by Hodson et al. [2, 3, 4, 5], 
Halstead et al. [6, 7], Qiu et at. [8], Sohn et ai [9] and 
Boyle et al. [10]. These investigations have helped 
identify /clarify the roles that factors such as the 
Reynolds number, free stream turbulence intensity, 
pressure gradient and curvature have in the genera- 
tion of losses. In particular, it has been determined 
that [1]: 

• At low to moderate Reynolds numbers there is 
a laminar region extending some distance from 
the leading edge. 

• The boundary layer may separate, particularly 
on the suction surface of the blade. Separation 
may occur in the form of a closed bubble, or 
as massive separation with no reattachment (re- 
sulting in large losses). The pressure surface may 
have cove separation, and small separation bub- 
bles may exist near the leading edge. The separa- 
tion bubbles often originate in transitional flow, 
while reattachment is usually in turbulent flow, 

• The interaction of incoming wakes with the 
boundary layer often creates a convected tran- 
sitional or turbulent patch, which is trailed by a 
“calmed” region. The calmed region is a relax- 
ation region between the patch and the laminar 
boundary layer. 

In parallel to the experimental investigations, there 
have been significant analytical efforts to improve the 
modeling of transition. Examples of such efforts in- 
clude the works of Mayle [11], Reshotko [12] and 
Gostelow et ai [13, 14]. These newer models show 
promise of providing accurate transition predictions 
over a wide range of flow conditions [15], although 
they have yet to be implemented into the numer- 
ical flow analyses used by the turbine design com- 
munity. Some recent computational investigations of 
interest include the works of Chernobrovkin and Lak- 
shminarayana [16], Kang and Lakshminarayana [17], 
Huang and Xiong [18] and Eulitz and Engel [19]. 

The focus of the current effort has been to use 
a viscous, unsteady quasi-three-dimensional Navier- 


Stokes analysis to study boundary layeT development 
in a two-stage low- pressure turbine. A two- layer alge- 
braic turbulence model, along with natural and bub- 
ble transition models, have been evaluated at both 
take-off and cruise operating conditions. The ge- 
ometry used in the simulations has been the sub- 
ject of extensive experiments [6, 7]. A recent inves- 
tigation by Dorney et ai compared the results of 
quasi-three-dimensional simulations with experimen- 
tal data for the same geometry [20]. The predicted 
time-averaged boundary layer quantities showed ex- 
cellent agreement with the experimental data, but 
the predicted unsteady results showed only fair agree- 
ment with the experimental data. It was hypothe- 
sized that unsteady effects associated with modified 
blade counts in the numerical simulations generated 
more unsteadiness than was observed in the exper- 
iments. The goal of the current work has been to 
quantify these effects of using more accurate blade 
counts, and to study the effects of post-processing 
techniques on the predicted results. 

ALGORITHM 

In the computational procedure the flow field is di- 
vided into two types of zones. O-type grids are used 
to resolve the flow field near the airfoils. The O-grids 
are overlaid on H-grids which are used to resolve the 
remainder of the flow field. The H-grids are allowed 
to slip relative to one another to simulate the relative 
motion between rotors and stators. The thin-layer or 
full Navier-Stokes equations are solved on both the 
O- and H-grids. The governing equations are cast in 
the strong conservation form. A fully implicit, finite- 
difference method is used to advance the solution of 
the governing equations in time. A Newton- Raph son 
subiteration scheme is used to reduce the lineariza- 
tion and factorization errors at each time step. The 
convective terms are evaluated using a third-order- 
accurate upwind-biased Roe scheme. The viscous 
terms are evaluated using second-order accurate cen- 
tral differences and the scheme is second-order accu- 
rate in time. Message Passing Interface (MPI) soft- 
ware has been implemented into the numerical anal- 
ysis to reduce the computation time for large-scale 
simulations. Details of the solution procedure and 
boundary conditions are discussed in Ref. [21] 

TURBULENCE MODEL 

The two- layer algebraic model based on the work 
of Baldwin and Lomax (BL) was used to model tur- 
bulence [22]. Several modifications were made to the 
original BL model based on previous experiences with 
compressor and turbine geometries: 
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ABS Model 


« The switchover location between the inner and 
outer models cannot move more than a specified 
number of grid points between adjacent stream- 
wise locations. This eliminates non-physical gra- 
dients in the turbulent viscosity near separation 
points. 


The ABS model determines the start of transition 
based on the momentum thickness Reynolds number 


Rc$ = 163 + exp 


\F(\e)- 


OM 

6.91 


Tu\ 


( 1 ) 


• A second derivative smoothing function is used 
on the turbulent viscosity field in separated flow 7 
regions. This also helps remove non-physical gra- 
dients in the turbulent viscosity in separation re- 
gions. 

• A cutoff value is imposed on the turbulent viscos- 
ity (nominally 1200 times the free stream lami- 
nar viscosity). 

The comparison of predicted and experimental in- 
tegral boundary layer quantities warrants discussion 
of the technique used to determine the location of the 
boundary layer edge in the simulations. The follow- 
ing steps, based on the work of Davis et ai [23], w r ere 
used to determine the edge of the boundary layer: 


where 

F{ As) = 6.91 + 12.75A, + 63.64(A e ) 2 A# < 0(2) 
F(A S ) = 6.91 + 2.48A« - 12.27(A,) 2 A* > 0 (3) 


This model is considered valid up to turbulence levels 
of Tn = 10%. The end of the transition region is 
calculated as 

Rel = 31.8(J?e*) 16 (5) 

In the region between the start and end of transition 
the intermittency function, a, is determined using the 
model developed by Dhawan and Narasimha. [25], 


1. determine the minimum value of |(7oo — u|, where 
Uqo is based on isentropic conditions 

2. correct the location based on where the local vor- 
ticity exceeds a specified limiting value 

3. within the new 7 range, determine where u is 
greater than Uoo 

The results presented in Ref. [20] for the current 
geometry exhibited larger variations in the unsteady 
displacement and momentum thickness distributions 
than w^ere observed in the experiments. One possible 
cause of the discrepancies w 7 as that the experimen- 
tal data acquisition system could measure bound- 
ary layer thicknesses to a maximum height of ap- 
proximately 1.2 mm, while there were no limitations 
placed on the boundary layer height in the simula- 
tions. Therefore, as part of post- processing in the 
current study the predicted boundary layer thick- 
nesses w 7 ere clipped at a maximum height of 1.2 mm. 


<j—l — exp (— 4.64£ 2 ) (6) 

Note, the intermittency function, wdiich has a value of 
<7 — 0 for laminar flow 7 and cr = 1 for turbulent flow 7 , is 
multiplied by value of the viscosity calculated in the 
turbulence model. The current implementation of the 
ABS model has been tested on both compressor and 
turbine geometries [26, 27]. 

Roberts* Model 

For cases involving separation bubbles the model 
developed by Roberts [28], and modified by Davis et 
ai [29], is used: 

Ree = 250001og 10 (1/ tanh(0.173205Tu)) (7) 

Instantaneous transition is assumed using the bubble 
model. 

GEOMETRY AND GRID 


TRANSITION MODELS 

The low Reynolds number environment in low- 
pressure turbines suggests that the flow may be tran- 
sitional. The natural transition model utilized in this 
study w 7 as developed by Abu-Ghannam and Shaw 
(ABS) [24], 


The test article used in this study has been studied 
extensively by Halstead et ai [6, 7]. The turbine is 
typical of those found in modern aircraft engines (see 
Fig. 1). The turbine design parameters are shown in 
Table 1. 

The experimental turbine contains 82 first-stage 
vanes, 72 first-stage rotors, 108 second-stage vanes 
and 72 second-stage rotors (i.e. f a 82-72-108-72 blade 
count). In the previous study a blade count of 78-78- 
104-78 was assumed, thus requiring the use of 3 first- 
stage vanes, 3 first-stage rotors, 4 second-stage vanes 
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RESULTS 


Quantity 

N1 

N2 

Rotor 

Solidity 

1.64 

1.68 

1.48 

Aspect Ratio 

3.97 

5.10 

3.87 

Chord, mm 

76.2 

59.7 

78.7 

Stagger, deg from axial 

39.8 

22.9 

34.9 

Camber, deg 

60.0 

102.4 

96.9 

No. of airfoils 

82 

108 

72 

Axial gap, mm 

25.4 

25.4 

- 


Table 1: Turbine design parameters. 


and 3 second-stage rotors in the simulations [20]. In 
the current study a 81-72-108-72 blade count ap- 
proximation has been used, requiring the use of 9 
first-stage vanes, 8 first-stage rotors, 12 second-stage 
vanes and 8 second-stage rotors in the simulations. 
For completeness, however, a set of 3-3-4- 3 simula- 
tions were computed using the boundary layer clip- 
ping technique described earlier. 

In keeping with the previous simulations, the 
O-grids for each blade passage contain 251 x 51 
(streamwisex tangential) grid points. The H-grids in 
the first-stage nozzle passage each contain 108 x 96 
grid points, the H-grids in the first-stage rotor pas- 
sage each contain 125 x 96 grid points, the H-grids in 
the second-stage nozzle passage each contain 120 x 96 
grid points and the H-grids in the second-stage rotor 
passage each contain 175 x 96 grid points. Thus, the 
complete grid topology contains 935,589 grid points 
(see Fig. 2, where every fifth grid point is shown for 
clarity). 

The average value of y + , the non-dimensional dis- 
tance of the first grid point above the surface, was 
approximately 1.0 for all four blade rows. The dimen- 
sions of the O-grid were arrived at based on the value 
of and the number of points within the boundary 
layer, while the dimensions of the H-grid were deter- 
mined by performing wake convection simulations in 
the absence of airfoils. 

The simulations were run on four processors of a 
Silicon Graphics Inc. (SGI) Origin200 workstation 
with 195-MHz processors. The average computa- 
tion time was 10.0 secs/iteration. Each simulation 
was run for 10 global periods, at 50,000 time steps 
per global period. A global period is defined as the 
second-stage rotor moving through a distance equal 
to 12 second-stage nozzle pitches or, similarly, the 
first-stage rotor moving through a distance equal to 
9 first-stage nozzle pitches. The 10 global periods 
enabled the efficiency, losses and integral boundary 
layer quantities to become time periodic. 


Two different operating points have been stud- 
ied: one at takeoff conditions (Point 5A) and one 
at cruise conditions (Point 5D). The parameters as- 
sociated with the two test conditions are outlined in 
Table 2. In accordance with the experiments, the free 
stream turbulence level (used in the transition model) 
was set at Tu — 3%. 



Test Point 

Quantity 

5A 

5D 

R c 

101,200 

34,300 

T n {°K) 

287.8 

287.8 

Pti(N/m 2 ) 

97,642 

100,839 

Pi (N/m~) 

96,514 

100,716 

Qii deg from axial 

0.0 

0.0 

fi(PPM) 

599.7 

185.1 

Ah 

0.1289 

0.0417 

Ps/Pt 1 

0.8614 

0.8614 


Table 2: Turbine design parameters. 


Point 5A 

Entropy contours are useful for tracking the con- 
vection of airfoil wakes. Figure 3 shows instanta- 
neous non-dimensional entropy contours for the 9- 
8-12-8 blade count. The stretching of the wakes be- 
tween the suction and pressure surfaces of adjacent 
airfoils is evident. In addition, it is possible to track 
the wakes from the first-stage nozzle downstream into 
the second-stage rotor passage. Instantaneous Mach 
contours are shown in Fig. 4. The passage-to-passage 
variations present in both Figs. 3 and 4 underscore 
the importance of using an accurate blade count ap- 
proximation. 

The predicted and experimental normalized blade 
loadings for the second nozzle and second rotor are 
shown in Figs. 5 and 6, respectively. The loading 
distributions were similar in both the 9-8-12-8 and 
3-3-4-3 simulations, therefore only the results from 
the 9-8-12-8 simulation are shown. The predicted re- 
sults exhibit good agreement with the experimental 
data, except in the suction surface leading edge region 
of the second nozzle. The differences in this region 
suggest the flow entering the second nozzle in the 
simulations has somewhat greater positive incidence. 
On the second rotor, both the numerical results and 
experimental data show an acceleration region near 
the suction surface leading edge followed by a zone 
of constant velocity, a second accelerat ion region and 
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finally a deceleration zone as the flow moves down- 
stream of the throat into the uncovered portion of the 
passage. 

Figure 7 contains the experimental and predicted 
(with clipped boundary layers) displacement thick- 
ness histories at 50% of the wetted distance on the 
suction surface of the second-stage nozzle over 4 ro- 
tor wake-passing periods. For comparison purposes, 
Fig. 8 displays the effects of the boundary layer clip- 
ping on the predicted 3-3-4- 3 results. The differ- 
ent blade count approximations have little effect on 
the predicted results, with both simulations showing 
the experimentally-observed increase in displacement 
thickness with the passing of the rotor wakes. The 
effects of clipping the boundary layer are more pro- 
nounced than the effects of blade count approxima- 
tion at this location (see Fig. 8). The clipping elim- 
inates the peaks observed in the results presented in 
Ref. [20], while preserving the character of the bound- 
ary layer in between wake passings. 

The experimental and predicted displacement 
thickness histories at 68% of the suction surface wet- 
ted distance are shown in Fig. 9, while the corre- 
sponding comparison of the clipped and undipped 
solutions are shown in Fig. 10. At this location, the 
distribution from the 9-8-12-8 simulation shows more 
unsteadiness than was predicted in the 3-3-4-3 sim- 
ulation, most notably just before the passing of the 
rotor wake. Both simulations indicate more rapid de- 
cay of the displacement thickness after the passing of 
the rotor wake than was observed in the experiments. 
Similar to the results at 50%), the effects of clipping 
are greater than the effects of the blade count ap- 
proximation (see Fig. 10). Note, while the unsteady 
results in Fig. 9 exhibit only fair agreement with the 
experimental data, the match is significantly closer 
than in the previous simulations without boundary 
layer clipping [20]. 

At 82% of the wetted distance on the suction 
surface the predicted displacement thickness history 
from the 9-8-12-8 simulation shows much closer agree- 
ment with the experimental data than the results of 
the 3-3-4-3 simulation (see Fig. 11). In particular, the 
9-8-12-8 simulation begins to mimic the gradual de- 
crease in the displacement thickness after the passing 
of the rotor wake. The effects of boundary layer clip- 
ping become more important approaching the trail- 
ing edge region (see Fig. 12). Similar to the results in 
Figs. 8 and 10 the clipping has only minimal effects 
on the histories in between the passings of the rotor 
wakes. 

Figure 13 illustrates the predicted and experimen- 
tal displacement thickness histories at 94% of the wet- 
ted distance on the suction surface. The results of 


both simulations exhibit only fair agreement with the 
experimental data, although the results of the 9-8- 
12-8 simulation lie slightly closer to the data. The 
discrepancies beteen the predicted results and the 
experimental data, however, are much smaller than 
those without boundary layer clipping (sec Fig. 14 
and Ref. [20]). The effects of the clipping appear 
to be two-fold. First the clipping makes the post- 
processing in the analysis more representative of the 
data acquisition in the experiments. Second, the clip- 
ping appears to account for the relaxation of the 
boundary layer, which is not included in the Baldwin- 
Lomax turbulence model. 

Comparisons of the experimental and predicted 
time- averaged suction surface displacement thickness 
distributions are shown in Fig. 15. Also included in 
Fig. 15 are the maximum and minimum displacement 
thickness distributions predicted in the simulations. 
The results of both simulations show close agreement 
with the experimental data. The 9-8-12-8 solution 
displays more unsteadiness, as indicated by the area 
between the minimum and maximum thickness dis- 
tributions. 

Point 5D 

The entropy and Mach contours at cruise condi- 
tions are similar to those shown in Ref. [20], and 
therefore are omitted here. 

The predicted and experimental displacement 
thickness histories at 50%< of the wetted distance 
along the suction surface of the second-stage nozzle 
are show in Fig. 16. Similar to the results at the 
takeoff conditions, the effects of blade count approx- 
imation are small at this location. The effects of the 
boundary layer clipping are also relatively small (see 
Fig. 17), although it must be noted that the unclipped 
distributions (in general) did not display the rapid, 
large amplitude variations observed at the cruise op- 
erating conditions [20]. The predicted results show 
only fair agreement with the data, as the decrease in 
the displacement thickness after the passing of the 
rotor wake is overestimated. 

Figure 18 contains the predicted and experimen- 
tal displacement thicknesses at 68% of the suction 
surface wetted distance. The results of the 9-8-12- 
8 simulation display better phase agreement with 
the data than the results of the 3-3-4-3 simulation. 
The results of both simulations show the double- 
peak nature of the distribution observed in the ex- 
periments. The effects of the boundary layer clip- 
ping at the 68% location, which are shown in Fig. 19, 
are greater than at the 50% location. The increas- 
ingly large effects of the clipping with wetted dis- 
tance is not surprising because the unsteadiness in 
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the boundary layer increases downstream of the tran- 
sition /separation points. 

At 82% of the wetted distance on the suction sur- 
face the predicted displacement thickness histories 
show fair agreement with the experimental data (see 
Fig. 20), although the results of the 3-3-4-3 simulation 
more closely resemble the shape of the experimental 
distribution. The clipping has a significant effect at 
the 82%) location where the solution displays two dis- 
tinct peaks and valleys per blade passing (similar to 
the experimental data), while the unclipped solution 
shows only one (see Fig. 21). 

Figure 22 illustrates the predicted and experimen- 
tal displacement thickness histories at 94% of the wet- 
ted distance on the suction surface. Although both 
simulations show good agreement with the experi- 
mental data, the shape of the distribution from 9- 
8-12-8 simulation more closely resembles that of the 
data. The effects of the boundary layer clipping are 
minimal at this location (see Fig. 23). 

Comparisons of the experimental and predicted 
time-averaged suction surface displacement thickness 
distributions are shown in Fig. 24. The results of both 
the 3-3-4-S and the 9-8-12-8 exhibit excellent agree- 
ment with the experimental data. The amplitude of 
the unsteady envelopes for the two simulations are 
similar, except near the leading and trailing edges. 
Near the leading edge the results of the 3-3-4-S simu- 
lation show more unsteadiness, while near the trailing 
edge the results of the 9-8-12-8 simulation show more 
unsteadiness. 

CONCLUSIONS 

A series of numerical simulations have been per- 
formed for flow through a two-stage low-pressure tur- 
bine. The simulations included two different operat- 
ing conditions (takeoff and cruise), two blade count 
approximations and two different methods for pro- 
cessing the boundary layer information. Comparisons 
between the predicted results and the experimental 
data indicate: 

• the effects of blade count, approximation increase 
with wetted distance along the airfoil suction 
surface (i.e., in regions where unsteadiness is the 
largest) 

• the method of obtaining/processing boundary 
layer information in the simulations must be con- 
sistent with the experimental techniques 

• the effects of boundary layer clipping were gen- 
erally more significant than the effects of blade 
count approximation 


• the predicted unsteady displacement thickness 
histories showed fair agreement with the experi- 
mental data, but much better agreement than in 
a previous study 

• the time-averaged displacement thickness and 
airfoil loading distributions show excellent agree- 
ment with the experimental data 

ACKNOWLEDGEMENT 

The work of the first author w T as supported by 
NASA Glenn Research Center under Grant NCC3- 
645. 

References 

[1] Simon, F. F. and Ashpis, D. E., “Progress in 
Modeling of Laminar to Turbulent Transition on 
Turbine Vanes and Blades,” International Con- 
ference on Turbulent Heat Transfer, 1996, also 
NASA TM 107180. 

[2] Hodson, H. P. } “The Development of Unsteady 
Boundary Layers on the Rotor of an Axial Flow 
Turbine,” AGARD CP-351, 1983. 

[3] Hodson, H. P., “Measurements of Wake- 
Generated Unsteadiness in the Rotor Passages 
of Axial Flow Turbines,” ASME Paper 84-GT- 
189, 1984. 

[4] Hodson, H. P., “Boundary Layer Transition and 
Separation Near the Leading Edge of a High- 
Speed Turbine Blade,” Journal of Engineering 
for Gas Turbines and Power , Vol. 107, 1985, pp. 
127-134. 

[5] Hodson, H. P., Huntsman, I. and Steele, A. B., 
“An Investigation of Boundary Layer Develop- 
ment in a Multistage LP Turbine,” Journal of 
Turbomachinery , Vol. 116, 1994, pp. 375-383. 

[6] Halstead, D. E., Wisler, D. C., Okiishi, T. H., 
Walker, G. J., Hodson, H. P. and Shin, H.-W., 
“Boundary Layer Development in Axial Com- 
pressors and Turbines - Part 1 of 4: Compos- 
ite Picture,” ASME Paper 95-GT-461, Houston, 
TX, 1995, also published in the ASME Journal 
of Turbomachinery. 

[7] Halstead, D. E., Wisler, D. C., Okiishi, T. H., 
Walker, G. J., Hodson, H. P. and Shin, H.-W., 
“Boundary Layer Development in Axial Com- 
pressors and Turbines - Part 3 of 4: LP Tur- 
bines,” ASME Paper 95-GT-463, Houston, TX, 
1995, also published in the ASME Journal of 
Tu rbom achtne ry. 


NAS A/TM— 2000-20991 1 


[8] Qiu, S. and Simon, T. W., “An Experimental 
Investigation of Transition as Applied to Low 
Pressure Turbine Suction Surface Flows,” ASME 
Paper 97-GT-455, Orlando, FI. 

[9] Solin, K. H., Shyne, R. J. and Dewitt, K. J., 
“Experimental Investigation of Boundarj r Layer 
Behavior in a Simulated Low Pressure Turbine,” 
ASME Paper 98-GT-34, International Gas Tur- 
bine and Aeroengine Congress and Exhibition, 
Stockholm, Sweden, June 2-5, 1998. 

[10] Boyle, R. J., Lucci, B. L., Verhoff, V. G., 
Camperchioli, W. P. and La, LL, "Aerodynam- 
ics of a Transitioning Turbine Stator Over a 
Range of Reynolds Numbers,” 98-GT-285, Inter- 
national Gas Turbine and Aeroengine Congress 
and Exhibition, Stockholm, Sweden, June 2-5, 
1998. 

[11] Mayle, R. E., “The Role of Laminar-Turbulent 
Transition in Gas Turbine Engines,” ASME Pa- 
per 91-GT-261, Orlando, FL, 1991. 

[12] Reshotko, E., “Boundary Layer Instability, 
Transition and Control,” AIAA Paper 94-0001, 
Reno, NV, 1994. 

[13] Gostelow, J, P., Blunder, A. R. and Walker, 
G. J., “Effects of Free-Stream Turbulence and 
Adverse Pressure Gradients on Boundary Layer 
Transition,” ASME Journal of Turbomachinery , 
VoL 116, 1994, pp. 392-404. 

[14] Gostelow, J. P., Melwani, N. and Walker, G. 
J., “Effects of Streamwise Pressure Gradient on 
Turbulent Spot Development, ASME Paper 95- 
GT-303, Houston, TX, 1995. 

[15] Boyle, R. J. and Simon, F. F., “Mach Number 
Effects on Turbine Blade Transition Length Pre- 
diction,” ASME Paper 98-GT-367, International 
Gas Turbine and Aeroengine Congress and Exhi- 
bition, Stockholm, Sweden, June 2-5, 1998, also 
accepted for publication in the ASME Journal 
of Turbomachinery. 

[16] Chernobrovkin, A. and Lakshminarayana, B., 
“Development and Validation of Navier-Stokes 
Procedure for Turbomachinery Unsteady Flow,” 
AIAA Paper 97-3281, Seattle, WA, 1997. 

[17] Kang, D. J. and Lakshminarayana, B., “Numer- 
ical Prediction of Unsteady Transitional Flow 
Due to Rotor Stator Interaction,” AIAA Paper 
97-2752, Seattle, WA, 1*997. 


[18] Huang, P. G. and Xiong, G., “Transition and 
Turbulence Modeling of Low Pressure Turbine 
Flows,” AIAA Paper 98-0339, Reno, NV, 1998. 

[19] Eulitz, F. and Engel, I\., “Numerical Investiga- 
tion of Wake Interact ion in a Low Pressure Tur- 
bine,” ASME Paper 98-GT-563, International 
Gas Turbine and Aeroengine Congress and Ex- 
hibition, Stockholm, Sweden, June 2-5, 1998. 

[20] Dorney, D. J., Ashpis, D. E., Halstead, D. E. 
and Wisler, D. C., “Study of Boundary Layer 
Devlopment in a Two-Stage Low-Pressure Tur- 
bine,” AIAA 99-0742, 37th Aerospace Sciences 
Meeting, Reno, NV, January 11-14, 1999. 

[21] Dorney, D. J., and Verdon, J. M., "Numerical 
Simulations of LInsteady Cascade Flow,” ASME 
Journal of Turbomachinery , VoL 116, No. 4, Oc- 
tober, 1994, pp. 665-675. 

[22] Baldwin, B. S., and Lomax, H., “Thin Layer 
Approximation and Algebraic Model for Sep- 
arated Turbulent Flow”, AIAA Paper 78-257, 
Huntsville, AL, January, 1978. 

[23] Davis, R. L., Hobbs, D. E., and Weingold, H. D., 
“Prediction of Compressor Cascade Performance 
Losing a Navier-Stokes Technique,” ASME Jour- 
nal of Turbomachinery , VoL 110, No. 4, October, 
1988, pp. 520-531. 

[24] Abu-Ghannam, B. J., and Shaw, R., “Natural 
Transition of Boundary Layers - The Effects of 
Turbulence, Pressure Gradient, and Flow His- 
tory,” IMechE Journal of Mechanical Engineer- 
ing Science , VoL 22, No. 5, 1980, pp. 213-228. 

[25] Dhawan, S., and Narasimha, R., “Some Proper- 
ties of Boundary Layer Flow During Transition 
from Laminar to Turbulent Motion,” Journal of 
Fluid Mechanics , VoL 3., 1958. 

[26] Dorney, D. J., “A Comparative Study of Four 
Algebraic Transition Models,” SAE Journal of 
Aerospace , VoL 103, Section 1, 1994, pp. 112- 
119. 

[27] Dorney, D. J., and Ashpis, D., “Study of Low 
Reynolds Number Effects on the Losses in Low- 
Pressure Turbine Blade Rows,” AIAA 98-3575, 
34th AIAA/ASME/SAE/ASEE Joint Propul- 
sion Conference, Cleveland, OH, July 13-15, 
1998, also to appear in the International Journal 
of Turbo and Jet Engines. 


NASA/TM — 2000-209911 


7 


[28] Roberts, W. B., “Calculation of Laminar Sepa- 
ration Bubbles and Their Effect on Airfoil Per- 
formance/ 1 AIAA Paper 79-0285, New Orleans, 
LA, January, 1979. 

[29] Davis, R. L., Carter, J. E., and Reshotko, E., 
‘"Analysis of Transitional Separation Bubbles on 
Infinite Swept Wings,” AIAA Journal , Vol. 25, 
No. 3, March, 1987, pp. 421-428. 

[30] Halstead, D. E., Wisler, D. C., Okiishi , T. H., 
Walker, G. J., Hodson, H. P. and Shin, H.-W., 
“Boundary Layer Development in Axial Com- 
pressors and Turbines - Part 4 of 4: Computa- 
tions and Analyses,” ASME Paper 95-GT-463, 
Houston, TX, 1995, also published in the ASME 
Journal of Turbomachine ry. 


NASA/TM— 2000-20991 1 


8 


Low Speed Research Turbine 


Airflow 



• Models low - pressure turbine 
blading of a commercial 
GEAE engine family. 

• Instrumentation 

- X hot - film probe in 

axial - tangential plane 

- midspan only 

- supports modified to 

mitigate prong vibration 

• Nominal Design Point: 

RPM=600 U mid =126 ft/s 

<)>=1 .03 \j/=1.17 

krotor=1-3 k n ozzle=0-75 

RessL =1 00k “ 500 k 


Figure 1: Experimental two-stage LPT turbine facility. 
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Figure 2: Computational grid for two-stage LPT turbine. 
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Figure 3: Instantaneous non-dimensional entropy contours - Point 5 A. 
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Figure 4: Instantaneous Mach contours - Point 5A. 
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Figure 7: Displacement thickness history, 50% s.s., 
Point 5A, nozzle- 2. 

Figure 5: Normalized loading on nozzle-2, Point 5A. 
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Figure 9: Displacement thickness history, 68% s.s. 3 Figure 11: Displacement, thickness history, 82% s.s. 3 
Point 5A, nozzle-2. Point 5A, nozzle-2. 



Figure 10: Displacement thickness history, 68%’ s.s., Figure 12: Displacement thickness history, 82% s.s., 
Point 5 A, nozzle-2. Point 5A, nozzle-2. 
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Figure 13: Displacement thickness history, 94% s.s., 
Point 5A, nozzle-2. 



Figure 15: Displacement thickness envelope, Point 
5A nozzle-2 



Figure 14: Displacement thickness history, 94% s.s., 
Point 5 A, nozzle-2. 



Figure 16: Displacement thickness history, 50%; s.s., 
Point 5D, nozzle-2. 
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Figure 21: Displacement thickness history, 82% s.s., Figure 23: Displacement thickness history, 94% s.s., 
Point 5D, nozzle-2. Point 5D, nozzle-2. 



Figure 22: Displacement thickness history, 94% s.s., Figure 24: Displacement thickness envelope, Point, 
Point. 5D, nozzle-2. 5 |) n0 zzIe-2 
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